PROCEEDINGS 

OF SCIENCE 



MKPH-T-10-29 
HIM-2010-04 



O 

(N 

> 
O 

in 



Systematic errors in extracting nucleon properties 
from lattice QCD 



in: 



> 
OO 

in 

m 



Stefano Capitani, Bastian Knippschild 

Institutfiir Kernphysik, University of Mainz, Becher Weg 45, 55099 Mainz, Germany 
E-mail: knippsch@kph . uni-mainz . de| 



Michele Delia Morte, Hartmut Wittig 

Institutfiir Kernphysik and Helmholtz Institute Mainz, University of Mainz, Becher Weg 45, 
55099 Mainz, Germany 



X 



Form factors of the nucleon have been extracted from experiment with high precision. However, 
lattice calculations have failed so far to reproduce the observed dependence of form factors on 
the momentum transfer. We have embarked on a program to thoroughly investigate systematic 
effects in lattice calculation of the required three-point correlation functions. Here we focus on 
the possible contamination from higher excited states and present a method which is designed to 
suppress them. Its effectiveness is tested for several baryonic matrix elements, different lattice 
sizes and pion masses. 
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1. Introduction 

The calculation of mesonic and baryonic matrix elements from lattice QCD has made good 
progress in the last few years. Simulations with fully dynamical quarks have reached the physical 
pion mass at large lattice sizes of L > 5 fm. For many mesonic matrix elements the overall uncer- 
tainties are of a few percent and no discrepancies between experiments and theory can be detected. 
Baryonic matrix elements have not yet reached this accuracy. For instance the axial charge or the 
electric form factor are not compatible with the experimental data, even for the largest volumes and 
the smallest pion masses [fi]|. 

In this study we will focus on the electro-magnetic and axial form factors of the nucleon. The 
matrix element of the vector current can be expressed by the Dirac- and Pauli form factors F\{q ) 
and F2{q 2 ) in the following way: 

(N(p\s , )\V,\N(p,s))=u(p\/)^F l (q 2 ) + i^^F 2 (q 2 )^u(p,s), (1.1) 

where \N(p,s)) is the nucleon ground state with momentum p and spin s, and u(p,s) is a Dirac 
spinor with mass m^. The momentum transfer is given by q = p — p' and a^ v = [y^, y v ]/2. The 
matrix element of the axial current can be expressed in terms of the axial form factor Ga(^ 2 ) and 
the induced-pseudoscalar form factor Gp{q 2 ): 

(N(p',s') \A^\N(p,s)) = u{p',s') ^y 5 G A (q 2 ) + y,^-G P {q 2 )^ u(p,s). (1.2) 

We will neglect Gp{q 2 ) here completely and focus on the electro-magnetic form factors and the 
axial form factor at zero momentum transfer, which corresponds to the axial charge g A - 
The observed discrepancy between the experimental and lattice data may be due to systematic 
effects. These can be lattice artifacts, finite volume effects, large pion masses and contaminations 
from excited states. We will focus on the excited state contributions here but plan to address the 
other effects later due to our lattice setup of various lattice sizes, lattices spacings and pion masses. 
Our ensembles for computing matrix elements were generated as part of the "Coordinated Lattice 
Simulations" (CLS) initiative For our measurements we use fully dynamical @(a) -improved, 
two flavor Wilson fermions where our solver is Schwarz preconditioned and deflation accelerated 
[^]. Here we will present data for one lattice spacing of 0.069 fm and two lattice sizes of 64 x 32 3 
and 96 x 48 3 . The pion masses vary between 300 and 900 MeV [f| g|. 

2. The standard plateau method 

The two-point function in Euclidean space time for a nucleon is defined as ^ : 

C 2 (p,t) = ^e- r P^ v (j v ( t ,x)f^0)) ^ |p~ V Tr [T p (-/// + m)} . (2.1) 

X P 

We use the interpolating field Jy = e ahc (u a Cy$d h )u C y to create a nucleon where C represents the 
charge conjugation matrix, greek letters Dirac indices and latin letters color indices. The energy 
E p of a nucleon with momentum p is related to its mass via the dispersion relation E 2 = m 2 + p 2 . 
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The factor Zg is the coupling strength of the baryon-state with the vacuum. For the polarization 
matrix F p we use \ ( 1 + 70 ) ( 1 — 1 ! /3 75 ) which projects the nucleon to positive parity and polarizes it 
in the z-direction. To optimize the overlap of the nucleon correlation function with the ground state 
Jacobi-smearing [Q] with HYP-smeared links in the spatial Laplacian [JsJ] is used. 
The computation of the three-point functions is more involved than for the two-point functions 
due to the necessity of computing extended propagators. For the vector and the axial vector form 
factors the following two diagrams contribute: 

My,O;t s ,p') = ^e s P'0 » —^ *,u > (2-2) 

L u (y,0-,t s ,p)=Z eSP '® o x • (2-3) 



The quark propagators are contracted at a fixed sink timeslice t s and the outgoing momentum p' is 
induced via a Fourier transformation. The initial momentum is always set to zero so the momentum 
transfer carried by the photon is q = —p'. This object is then used as a source for a new set of 
inversions to create a so called extended propagator [§J. 
The three-point function can be constructed from the extended propagators as 

C 3 (q,t,t s ) = 1 £ i Tr[r p (Z u (0,y)±'L d (0,y))O(y)S(y,0)] e®, (2.4) 

y 

where the plus sign corresponds to the isoscalar form factor and the minus sign to the isovector 
form factor. Here we focus on isovector form factors where no disconnected contributions arise. 



The usual quark propagator S(y, 0) closes the Feynman diagrams of eq. ^2| and [1J and the operator 
O is inserted at timeslice t. 

At the hadronic level and for large Euclidean time the correlation function can be written as 

C 3 (q,t,t s ) = 2>- m *W^) Z | /^ Tr [T p u{p\s') (N(p', S f ) \O\N(0,s))u(0,s)] , (2.5) 
s,s' V E P 

where the matrix element (N(p',s ! ) \O\N(0,s)) can be expanded in terms of the form factors as 
shown in eq. [O] and |1.2[ We use in our simulations the local currents Vu(x) = x y(x)y fl x i , (x), 
= ^ '( x )757/i*^ '(x) where ^(x) is a u- or J-quark spinor. The local currents need to be renor- 
malised and in the case of the electro-magnetic form factors this can be done by requiring Ge (0) = 1 
which imposes charge conservation. For the axial current we use the non pertubatively computed 
renormalisation constant Za from JTo|]. 

Matrix elements computed with & (a) -improved Wilson fermions are not automatically &{a)- 
improved. An improvement term has to be added to the axial current, but we are interested in 
the axial charge only which is extracted from a spatial component at zero momentum of the axial 
current and so the forward matrix element of the improvement term vanishes. Here we neglect the 
improvement term for the vector current but it will be included in a later stage of our analysis. 



The standard way to extract form factors from Euclidean correlation functions like eq. 2.5 is to use 
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ratios between two- and three-point functions. Our particular choice is [ffllj: 



R(q,t,t s ) 



C 3 (q, t , t s ) / C 2 (q, t s - t)C 2 (0, t)C 2 (0, t s ) 



(2.6) 



C 2 (0,t s ) V C 2 (0,t s -t)C 2 (q,t)C 2 {q,t s ) 
This ratio cancels the exponential factors and gives usually long plateaux with small errors. With 
our choice of the polarization matrix we get access to the following quantities: 

H'2 



Z v -Re[R(q,t,t 



lm + E q j ^ 



2£„ 



ZvRe^f.f 



2E q (E q + m 



2m(m + Epi) 
{F,{q 2 )+F 2 {q 2 )} 



Z A -\m[R(q,t,t s ) 



75 75 J 



Im+E, 

^En 



• ! >G A (q 2 )- 2 -^G p (q 2 : 



in 



(2.7) 
(2.8) 
(2.9) 



The masses and energies in the pre-factors of the bare ratios can be determined from the two-point 
function of the nucleon. In principle the form factors can also be extracted from other components 
of the ratio but these turned out to be too noisy within the standard approach. 
The plateau method should give the ground state value of the matrix elements, but higher state con- 
tributions can lead to wrong plateau values. In the left panel of fig. [I] we show the connected part of 
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Figure 1: Examples for the standard method for Vo and gA 

the zero component of the isoscalar vector form factor for three different source/sink combinations 
on a 64 x 32 3 lattice with a pion mass of m K = 550 MeV. The black points are from a point source, 
the purple ones from a Jacobi-smeared source and the red ones from a smeared source/sink com- 
bination. All three graphs give reasonable plateaux yet they are not compatible within statistical 
errors. 

In the right panel the ratio for the axial charge for three different sink positions t s and a Jacobi- 
smeared source on a 64 x 32 3 lattice with a pion mass of m n = 415 MeV is shown. The ratio 
should be independent of the position of the sink but it is obviously not. Finally we cannot con- 
clude that the real plateau value is obtained even for the largest sink timeslice t s = 15. 

3. The summation method 

We will now argue that the large dependence of the ratios on t s and different smearing levels 
is mainly induced by higher order corrections and propose a method to reduce them. The ratio for 
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an arbitrary operator can be written as: 

R{q,t,t s ) =R G +0(e- At ) + (e- A '( f <- f >) . (3.1) 

It can be expressed in terms of the value of the ground state Rq and the exponentially suppressed 
excited state contributions. The energy gaps A and A' can be different for different quantum num- 
bers at source and sink. Assuming these gaps come from a two-pion state these contributions are 
large for small pion masses and small t s . 

To deal with these effects an alternative method can be used known as the summation method [|l2|] . 
The basic idea is to sum the ratio in t up to t s . Doing so results in 

t R{q,t,t s ) = R G ■ t s + c(A,A') + (t s e-^) + Le^A , (3.2) 

where again only the first excited state is taken into account. The ground state can be extracted 
from the slope of a linear function in t s . Also here an exponantially suppressed higher order cor- 
rection survives but in contrast to the standard method it is only ? s -dependent and so much smaller. 
Fig. H illustrates how well the summation method works. We plot the zero component of the con- 
nected part of the isoscalar vector current for four different momentum transfers. The behavior is 




Figure 2: The summation method for Vb and different momenta 



linear as expected and no deviation from the linearity can be detected. This means that even for the 
smallest t s the excited state contributions are depleted. 

In order to compare the summation method directly with the standard method we show in fig. || 



1.8 
1.6 
1.4 
1.2 
1 

0.8 
(16 
0.4 
0.2 



: ^ Tt {{{ t t' 






v I 

■ 






1 * ¥ ^ 1 ' PP ^ 
PS,N=150 □ 
SS,N=150 v . 
summation method for PS 



10 12 14 16 18 



E 1.4 



summation method 



t s = 10 
t = 15 



***** i 



* A ^ a 



2 4 6 8 10 12 14 16 
t/a 



Figure 3: Examples for the standard method and the summation method for Vq and g& 
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the graphs of fig. [I] again, but now including the values from the summation method. In the left 
panel the purple line represents the summation method using a correlation function with a smeared 
source. The thick line in the middle corresponds to the mean value and the thinner outer lines show 
the statistical error. The summation method depletes the higher state contribution in both cases at 
the expense of a larger statistical error. In the right panel the ambiguity in t s is completely resolved 
by using the summation method (purple lines) and the trend seems to end there. 



4. Results 

Our results for the axial charge and the vector form factors are preliminary. As mentioned 
before the improvement term for the vector current is not implemented yet and the scale is not set 
to its final value so the momenta and pion masses may change. All values were extracted with 
the summation method. In fig. ^| we show our results for the Dirac and the Pauli form factors. 
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Figure 4: Preliminary results for the Dirac and the Pauli form factors extracted with the summation method. 

While our values for the Pauli form factor for the largest lattice with the smallest pion mass is in 
good agreement with the experimental curve the Dirac form factor shows a different slope. So the 
excited state contributions seem not to be the only effect which causes the discrepancy between 
lattice results and experiment. 

In fig. |5| we show our data for the axial charge for different lattice sizes and pion masses. Data from 



1.4 
1.3 
1.2 
1.1 
1 

0.9 
0.8 
0.7 



PRELIMINARY 



V=(1.7fmK 
V=(2.2fm), 
V=(3.3fm) 3 
RBC/UKQCD 
experiment 



■V ( GeV ) 

Figure 5: Preliminary results for the axial charge extracted with the summation method. 



the RBC/UKQCD collaboration [13] are shown for comparison. In contrast to other collaborations 
we do not see a downward trend for small pion masses. The impact of the summation method is 
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larger for small pion masses because the energy gaps A and A' decrease and the influence of excited 
state contributions increases. 

5. Summary and outlook 

We showed that the control of excited states is crucial and although plateaux of baryonic ma- 
trix elements extracted with the standard method may look reasonable this method is insufficient to 
ensure the absence of excited state contributions. Smearing appears not to be enough to suppress 
these contributions so we introduced an alternative method, the summation method. It is promis- 
ing to help with excited state contributions but more inversions are needed and the statistical error 
grows compared to the standard method. To minimize these drawbacks this method will be opti- 
mized and tuned. 

To be able to extrapolate the matrix elements to the physical point simulations with smaller pion 
masses and different lattice spacings (/3 = 5.2,5.5) are being analysed. In the analysis of the vector 
form factors the G(a) improvement terms will be included, and finally we are also interested in the 
full axial form factor Ga(^ 2 ) and the induced pseudo scalar form factor Gp(q 2 ). 
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